Accepted to ApJ on 201 1 September 5 

Preprint typeset using KTgX style emulateapj v. 1 1/10/09 



O 

(N 

Cu- 
te 

in 
in 



X 

6 
to 



> 
oo 

o 
o 



X 



NEUTRON STAR STIFF EQUATION OF STATE DERIVED FROM COOLING PHASES OF THE X-RAY BURSTER 

4U 1724-307 

VALERY SULEIMANOV, 1,2 JURI POUTANEN, 3 MIKHAIL REVNIVTSEV, 4 AND KLAUS WERNER 1 
1 Institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University, Sand 1, 72076 Tubingen, Germany; 

[suleimanov,werner] @ astro.uni-tuebingen.de 
2 Kazan Federal University, Kremlevskaja str., 18, Kazan 420008, Russia 
3 Astronomy Division, Department of Physics, PO Box 3000, FI-90014 University of Oulu, Finland; juri.poutanen@oulu.fi 
4 Space Research Institute (IKI), Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia 

Accepted to ApJ on 2011 September 5 

ABSTRACT 

Thermal emission during X-ray bursts is a powerful tool to determine neutron star masses and radii, if the 
Eddington flux and the apparent radius in the cooling tail can be measured accurately, and distances to the 
sources are known. We propose here an improved method of determining the basic stellar parameters using 
the data from the cooling phase of photospheric radius expansion bursts covering a large range of luminosities. 
Because at that phase the blackbody apparent radius depends only on the spectral hardening factor (color- 
correction), we suggest to fit the theoretical dependences of the color-correction versus flux in Eddington units 
to the observed variations of the inverse square root of the apparent blackbody radius with the flux. For that 
we use a large set of atmosphere models for burst luminosities varying by three orders of magnitude and for 
various chemical compositions and surface gravities. We show that spectral variations observed during a long 
photospheric radius expansion burst from 4U 1724-307 are entirely consistent with the theoretical expectations 
for the passively cooling neutron star atmospheres. Our method allows us to determine both the Eddington flux 
(which is found to be smaller than the touchdown flux by 15%) and the ratio of the stellar apparent radius to the 
distance much more reliably. We then find a lower limit on the neutron star radius of 14 km for masses below 
2.2M©, independently of the chemical composition. These results suggest that the matter inside neutron stars 
is characterized by a stiff equation of state. We also find evidences in favour of hydrogen rich accreting matter 
and obtain an upper limit to the distance of 7 kpc. We finally show that the apparent blackbody emitting area in 
the cooling tails of the short bursts from 4U 1724-307 is two times smaller than that for the long burst and their 
evolution does not follow the theory. This makes their usage for determination of the neutron star parameters 
questionable and casts serious doubts on the results of previous works that used for the analysis similar bursts 
from other sources. 

Subject headings: radiative transfer - stars: neutron - X-rays: bursts - X-rays: individual (4U 1724—307) - 
X-rays: stars 



1. INTRODUCTION 

Studies of the thermal emission from neutron stars 
(NSs) have been u sed extensively to determine their 
masses and radii (e.g. IDamen et alJll990t Ivan Paradiis et alJ 
19901: ILewin et all 119931: iRutledge et alj|2002t iHeinke etai] 
20061 TWebb & Barret! 120071: lOzel et all 120091: iGuver et all 



2010a]), which can provide constraints on the properties of 



the matter at supranucle ar densities (Hae nsel et al.1 12007b 
Lattimer & Prakash 2007). Thermonuclear X-ray bursts at 
NS surfaces and, particularly, the photospheric radius expan- 
sion (PRE) bursts exceeding t he Eddington limit at high-flux 
phases ( Kuulk ers et al.1 120031) are excellent laboratories for 
such detailed studies. The observed Eddington flux gives 
a constraint on the NS mass-radius relation. This method, 
however, suffers from the uncertainty in determination of 
the exact moment when the luminosity reaches the Edding- 
ton limit at the NS surface, which is often assumed to coin- 
cide with the moment of "touchdown", when the measured 
col or temperature is highest and the apparent radius is small- 
est (IDamen et al.lll990f) . This interpretation is uncertain, as 
the t ouchdown flux typical ly coincides with the maximum 
flux dGalloway et al.ll2008bl) . while the latter is expected to be 
larger than the surface E ddington flux by the redshift factor 
l + z (ILewin et al.lll993l) . 
The second constraint can be obtained from the apparent ra- 



dius of the NS at late stages of the burst. This method also has 
systematic uncertainties related to the color-correction factor 
f c — T c /T e s (the ratio of the color temperature to the effec- 
tive temperature of the star), which is a function of the burst 
luminosity. 

To minimize the theoretical uncertainties in modeling the 
burst a tmospheres at nearly Ed dington luminosities, we pro- 
posed ( Suleima nov et al.ll"201 ll) to use the whole cooling track 
and check that the evolution of the blackbody normalization 
with flux is consistent with the theoretically predicted evolu- 
tion for a passively cooling neutron star. 

The unique determination of mass and radius requires one 
more constraint. Knowing the distance to a burster breaks the 
degeneracy, and therefore bursters located in globular clusters 
would be of interest. For the analysis in the present paper 
we choose 4U 1724-307 which resides in the globular clus- 
ter Terzan 2. We analyze three PRE bursts from that source 
and show that they follow different cooling tracks. We also 
speculate on the reasons for such a discrepancy. We demon- 
strate that the apparent blackbody radius for the long PRE 
burst evolves according to the predictions for the passively 
cooling neutron star with a constant apparent surface and use 
these data to determine the NS parameters. 
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Figure 1. Spectra of the persistent emission before the long burst (open 
squares) and one of the short bursts (filled circles). 



2. DATA 



2.1. 



X-ray bursts from 4U 1724-307 

For our analysis we have used the data from the PCA 
spectrometer of Rossi X-ray Timing Explorer (RXTE), be- 
cause it provides us the maximum possible number of pho- 
tons within the small time intervals, during which the spec- 
trum of the X-ray burst does not vary much. RXTE/PCA data 
were analyzed with the help of HEASOFT package (version 
6.6.1). Response matrices were generated using task PCARSP 
(v.10.1) of this package. The background of PCA detectors 
was estimated with the help of CM_bright_VLE model. 
In order to account for the uncertainties reflecting the qual- 
ity of the RXTE/PC A response matrix, a 1 % systematic error 
dJahoda et alj |2006) was added in quadrature to the statistical 
error in each PCA energy channel. Th e spectral fittin g was 
performed using the XSPEC 1 1 package dArnaudll 19961) . 

RXTE observed thr ee PRE bursts from 4U 1724-307 
(lGallowayetalJl2008al) . A lon g (>150 s) PRE bu rst was 
recorded on November 8, 1996 ( Mol kov etaTl 12000). Two 
short PRE bursts were observed on Feb 23 and May 22, 2004. 
The quiescent emission around these three bursts was signif- 
icantly different. The long burst happened when the source 
was in the so called hard/low state with a luminosity of a few 
percent of the Eddington luminosity and the X-ray emission 
was formed in an optically thin medium. The short bursts 
were observed during the high/soft state, with the X-ray emis- 
sion characterized by the optically thick accretion disk and the 
boundary/spreading layer from the neutron star surface. The 
quiescent spectra of 4U 1724-307 are presented in Fig.Q] All 
spectra presented here and below were fitted assuming an ab- 
sorption column density of Nh = 10 22 cm -2 corresponding 
to the best-fit value for the persistent spectrum. But fitting 
using Nh as a free parameter was also performed. 

X-ray bursts often have thermally looking spectra which 
can be approximated by a blackbody (Gallo way et alj |2008a). 
The fits of the burst spectra with an absorbed blackbody 
model give, however, rather high reduced x 2 /d.o.f^2-3. This 
is not very surprising because the photon statistics is high 
with the errors being dominated by the systematics. Theo- 
retical models of the NS atmosphere during X-ray bursts do 
not predict exactly blackbody spectral shapes, but the resid- 
uals in the 3-20 keV RXTE e nergy band are expected at the 
level of 4-5% (see Fig. 8 in Suleima nov et al.H2.Ql 11) . The 
residuals between the blackbody model and the data of X-ray 
bursts of 4U 1724-307 have similar amplitudes (see Figs. [5] 
and [3). Adopting the systematic uncertainty of the model 
at the level of 3-5% allows us to obtain acceptable x 2 val- 
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Figure 2. Evolution of the long burst spectra at early decline burst phases 
with the corresponding best-fit blackbody models (solid curves). The corre- 
sponding time points are shown by arrows in the middle panel of Fig.f4] 
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Figure 3. Same as Fig. ff] but for a short burst. 



ues for all spectra at the cooling tails of the analyzed X-ray 
bursts, while simple blackbody analytical model allows us to 
correctly represent the general shape of the spectrum in 3-20 
keV energy band. Therefore in our subsequent analysis we 
use simple blackbody models for approximating the spectral 
shapes and compare the obtained parameters with those, ob- 
tained by fitting the spectra produced by the full theoretical 
NS atmosphere model in the same energy band. 

The best-fit parameters of the employed blackbody model 
are the color temperature Tbb and the normalization constant 
K = (i?bb[km]/£>io) 2 (the blackbody radius is measured in 
km and the distance D to the source in units of 10 kpc). These 
can be combined to estimate the observed bolometric flux F . 

The evolution of the fitted parameters during the bursts are 
shown in Fig.|4] The time is normalized to the individual 
flux decay timescales r and shifted to allow an easy compar- 
ison between the bursts. The time point, when the flux and 
the color temperature Tbb reach their maximum values and 
the normalization K has a minimum, is usually named the 
"touchdown" point (shown by the arrow in the upper and mid- 
dle panels of Fig.|4]for the long burst). The maximal fluxes of 
the short bursts are appreciably smaller than that for the long 
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Figure 4. Evolution of the observed blackbody fluxes, color temperatures 
and normalizations K = R+/Df for three bursts from 4U 1724-307 in 
November 8, 1996 (black circles), February 23, 2004 (blue diamonds) and 
May 22, 2004 (red triangles). The time variable is normalized to the char- 
acteristic decay time r which is equal to 18.5, 2.54 and 3.52 s for these 
three bursts, respectively. For the long burst, the zero time corresponds to 
the touchdown point (marked by an arrow in the upper and middle panels), 
while for the short burst the light curves were shifted so that the coolingtails 
coincide. The arrows show the times when the spectra shown in Figs. |2]and 
fl]were collected. 

burst, and their normalizations do not show a significant rise 
at the early stages of the bursts indicating that the NS photo- 
sphere has not substantially expanded. 

The most serious difference between the short and the long 
bursts is the normalization value at the flux decay phase, the 
apparent area in the long burst in approximately two times 
larger (Pig.|5j. The natural explanation for that is the presence 
of the optically thick accretion disk which in the soft state 
blocks a significant part of the NS apparent surface, while the 
hot, optically thin, transparent accretion flow in the hard state 
does not affect much the NS apparent area. In addition to that 
effect, there could be additional differences in the physical 
conditions of the NS atmosphere (and thus in spectral harden- 
ing factors), as the boundary/spreading layer in the soft state 
forces the NS atmosphere to r otate close to Keplerian velocit y 
(Inoga mov & Sunvaevlll999t lSuleimanov & Poutanen 2006). 
We note that similar differences between short and long bursts 
were a l so found in anothe r X-ray bursting NS (Gallowav et al. 
l27MilZnan7et^l2lH0h . 

2.2. Distance to 4U 1724-307 

4U 1724-307 resides in the globular cluster Terzan 2. 
The distance t o that J=7.5± 0.7 kpc was measur ed by 
iKuchinski et"all (119951) . while lOrtolani et all (119971) give 
D=5.3±0.6 kpc, if Ry = A V /E(B - V) = 3 6 (more suit- 
able for red stars, see iGrebel & Roberti fl995l) . or 7.7±0.6 
kpc if Ry = 3.1. To cover all possibilities, we assume further 
a flat distribution from 5.3 to 7.7 kpc with Gaussian tails of 
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Figure 5. Evolution of the bursts in the flux - temperature plane for the three 
bursts from Fig.f4] The curves of constant describing the decay phase 
are also shown. The blackbody apparent area for the short bursts (dotted 
curve) is a factor of two smaller than the corresponding area describing the 
long burst (dashed curve). 

lcr=0.6 kpc on both ends. 

3. METHOD 

3.1. Basic relations 

Here we briefly present some well known relations between 
observed and real physical NS parameters, which arise due 
to the gravitational redshift and the light bending, together 
with some other important equations. The observed luminos- 
ity Lqo, effective temperature Too and apparent NS radius i?oo 
are connected with the luminosity at the NS surface L, the ef- 
fective temperature measured at the surface T c g, and the NS 
circumfere ntial radius R and mass M by the following rela- 
tions dLewin et alJI 19931) 

Loo = - L Too = -^1 R oa = R{l + z), (1) 



(l + z) 2 
with the redshift factor 

l + z = 



l + z' 



'1 - 2GM/Rc 2 



-1/2 



(2) 



The gravity g on the NS surface is larger in comparison with 
the Newtonian case due to the general relativity effects 

GM . 



9 = 



R 2 



■(l + z), 



therefore, the Eddington luminosity is larger too: 

AnGMc , , , 



Edd- 



(3) 



(4) 



Here TEdd is the maximum possible effective temperature 
on the NS surface, k c — 0.2(1 + X) cm 2 g _1 is the elec- 
tron (Thomson) scattering opacity, and X is hydrogen mass 
fraction. The observed Eddington luminosity is smaller for 
higher z inGMc 1 

-kEdd.oo = 



K c 1 + Z 

This is related to the observed Eddington flux 

„ _ LEdd,oo _ GMc 1 
FEdd ~ l^D 2 " 
and the Eddington temperature 



, D 2 l + z 



(5) 



(6) 
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1/4 



-Roo 

IT 



-1/2 



(7) 
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which is the effective temperature corresponding to the Ed- 
dington flux at the NS surface corrected for the gravita- 
tional redshift. We note here that the e lectron scattering 
opacity decreases with temperature (e.g. iPaczynskH 119831: 
iPozdniakov et alJl 19831) and at the typical temperature of ~3 
keV in the upper layers of the luminous NS atmospheres is 
reduced by about 7%. This affects the luminosity where the 
Eddington limit actually is reached. 

3.2. Atmosphere models and color correction 

Numerous computations of X-ray bursting NS at- 
mosph er es (lLondon et all 11984 1 1986b lEbisuzakil 11987 



Made il [199 lb iPavlov et all 1 1 99 lb IZavlin et all 
Suleimanov & Pout anenl 12006b ISuleimanov et al.1 1201 1 
show that their emergent spectra at high luminosities are 
close to diluted blackbody spectra due to strong energy 
exchange between high energy photons and relatively cold 
electrons at NS surface layers (Compton down-scattering) 

F E « w B E {T C = f c T cS ), (8) 

where f c is the color correction (or hardness) factor and w is 
the di lution factor, which at hi gh luminosities is very close to 
II f c 4 (ISuleimanov et al.ll20TTb . 

Spectra observed from the X-ray bursting NSs are close to 
thermal and usually they are fitted by a blackbody with two 
parameters: the observed color temperature Tbb and the nor- 
malization K = (i?bb(km)/Z?io) 2 . It is easy to find the rela- 
tions between various temperatures 

bb — Jc -Loo — ./ C 



Tec — fc 1 — 1 

1 + z 1 + z 
The observed blackbody flux is then 

£?2 r>2 



(9) 



F = <Ts B T b 4 b -Jf = <J SB T^ (10) 

and we can find the relation between the normalization and 
the NS radius 

2 n i ^2 D2 i 

(11) 

ft D 2 ft 



R bb 



R 2 (1 + z) 2 



D 2 D 2 

These formulae can be transformed to the relation b e tween 
color correction and normalization dPenninx et al.l 119891: 
Ivan Paradiisetal .1119901) : 



K 



-1/4 



f c A, A = (R 00 [km]/D 



10) 



-1/2 



(12) 



A combination of A and i^Edd gives the Eddington tempera- 
ture: 



Edd.oo 



= 1.14 X 10 8 AR 



1/4 



_ 7 K = 9.8lA^/ 4 dj . 



Edd,-7 iv — a ' ui r Edd,-7 ^V, 

( 13 ) 

where F Ed d -7 = ^Edd/10 7 erg cm 2 s 1 . 

A detailed comparison of the theoretical models with the 
data requires the knowledge of the run of the color cor- 
rection with flux. Previous models covered the range of 
lumin osities very sparsely. Using our recently developed 
code (Sulei manov & Poutanenll2~006t ISuleimanov et al.ll2006l: 
Suleimanov & Wernerl l2007l) . we have computed a very de- 
tailed set of model s with the luminosity va rying by three or- 
ders of magnitude (Suleimanov et al. 20Tlh. 

An atmosphere model is fully defined by the surface gravity 
g, chemical composition, and the ratio of the luminosity to the 
Eddington luminosity I = L/Lec\c\- The last parameter also 
relates the effective temperature to the Eddington temperature 
at the NS surface: 




T, 
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l^T, 



Edd- 



(14) 
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Figure 6. Color- correction factor as a function of the NS luminosity 
ISuleim anov et al. 2011). The curves correspond to atmospheres of differ- 
ent chemical composition: pure hydrogen (red), pure helium (pink); the 
blue curves are for models with solar H/He composition plus solar abun- 
dance of metals Z = Zq (dotted blue) and subsolar metals Z = 0.3^0 
(solid blue). The surface gravity is taken to be g = 10 14 cm s — 2 . The 
dashed curve shows the results for a hydrogen atmosphere at larger gravity of 
log 9 = 14.3. 



We considered various chemical compositions (pure hydro- 
gen, pure helium, solar mixture of hydrogen and helium with 
various metal abundances) and three surface gravities log g = 
14, 14.3, and 14.6. The (redshifted) radiation spectrum from 
the NS atmosphere was then fitted with a diluted Planck func- 
tion in the 3-20 keV energy band (i.e. the range observed by 
RXTE) to determine the color-correction factor / c (see Fig. 
[6]). The behaviour of / c at relatively high I depends mainly 
on the hydrogen abundance X and very little on the surface 
gravity and metal abundance. 

3.3. Determining M and R using the touchdown method 

In an ideal situation if the observed X-ray emission comes 
indeed from the passively cooling fully visible NS and the dis- 
tance to the source is known, we can determine NS mass M 
and radius R from two observables: the Eddington flux given 
by Equation (O and the NS apparent blackbody size in the 
coolin g tail, or quantity A = K^ 1 / 4 / f c (see e.g. lLewin et alj 
[19931) . The latter is related to the apparent size of the NS 
through the color correction, Roo = f 2 Rbb, and / c is as- 
sumed to be known in the burst tail from the theoretical con- 
siderations. While this method was proposed long time ago, 
only recently strong claims appeared in the literature that it 
actually can be used for determining accurately paramet ers of 
three bursting NSs dOzel et alJl2009T:lGuver etal]|2010allbh . 

Using the approach advertised in the aforementioned pa- 
pers, one has to determine the Eddington flux from observa- 
tions. For PRE burst s it was assumed t hat it is reached at the 
"touchdown" point (Damen et al. 1990), when the color tem- 
perature is highest and the apparent blackbody area lowest. 
The color correction factor f c at the late c ooling phases of 
the PRE bursts was taken to be close to 1 .4 (lOzel et al.ll2009b 
Giiveretal. 2010a b ) based of the models by iMadej et alJ 
(2004) and Majczyna et al. (2005) (see discussion below). 
The observables can be then transferred to the constraints on 
M and R (see Fig. [7). We will further call this approach the 
"touchdown method". 

From the Eddington flux estimate we have (see dotted 
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Figure 7. Constraints on M and R from various observed values. The solid 
curve giv es th e relation obtained from the Eddington temperature given by 
Equation )13t . thick dotted curve is for the Eddington flux given by Equa- 
tion |6}, thick dashed curve is for j4=const. If the assumed distance is too 
large, there are no solutions (the corresponding curves for F Edd =const and 
A=const shown by thin lines do not cross). 
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D 2n c D 2 F Edd _! 
R = 5 u 1 



1-1/2 



14.138 km (1 + X) Dj F Edd ,_ 7 u' 1 (1 



and the mass is found using 

M 



R 



M Q 2.95 km 



(16) 

(1 + z)~ 2 and 



where the compactness u = Rs/R = 1 
i?s = 2GM j c 2 is the Schwarzschild radius of the NS. A mea 
surement of A gives another constraint: 



R = i?ooVT 



= fl DVK = D 10 A- 2 VT^Ukm. 

(17) 

Combining with the parametric expression (TTBT l for the mass, 
we get the second relation between M and R shown by the 
dashed curves in Fig. [7] 

The Eddington temperature given by Equation ( fT3l is in- 
dependent of the uncertain distance to the source and can be 
used to express the NS radius through the observables and u: 



R 



2KeC7SBT EdcLc 



u(l-u) 3 / 2 , 



(18) 



and the mass is then found via Equation ([Tol l. The correspond- 
ing relation between M and R is shown by the solid curve in 

Fig.m 

All three curves cross in one or two points (see Fig. 13 if the 
quadratic equation 

«(l-u) = 2K lf„f E f = 14.138 (l+X)D 10 F Edd ,_ 7 A 2 , 



'Kfl 



(19) 



which follows from Equations ( fT3T > and ( [T7| >, has a real so- 
lution for u (see e.g. ISteiner et all l2010h . This happens if 
u(l — u) < 1/4 and the distance then should satisfy the fol- 
lowing inequality 



Kf 2 c 3 



0.177 



8 i*Edd (1 + X) A 2 F Edd ,- 7 



kpc. 



(20) 

In the opposite case, there is no physical solution for M and 
R for given observables. 



As we mentioned above, this method of determination of 
M and R works in an ideal situation. There are a few prob- 
lems with this approach. First, the relation of the Eddington 
flux to the touchdown flux is not clear. The reduction of the 
electron scattering opacity at high temperatures increases the 
true Eddington limit by about 7% above that given by Equa- 
tions (0 and ©. Also if we believe that the X-ray burst lumi- 
nosity is equal to the Eddington luminosity during the expan- 
sion phase of the PRE burst, the observed luminosity has to 
decrease when t he photospheric ra dius decreases according 
to Equation (0 dLewin et alj|1993l) . In reality the observed 
luminosity in 4U 1724-307 increases when the photospheric 
radius decreases (see Fig. |4)l. This implies that the ratio of 
the luminosity to the Eddington limit at the photosphere has 
to increase with decreasing photospheric radius (i.e. increas- 
ing redshift z). A combination of this dependence with the 
gravitational redshift effect predicts then that the observed lu- 
minosity reaches the maximum when the photospheric radius 
is larger than R and that maximum is larger than expected for 
the Eddington luminosity at the surface. 

Second, the assumption of f c s» 1.4 in the cooling tail 
is very uncertain. Thi s assumption is based on X-ray burst 
atmosphere models of iMajczyna et all (120051) . who claimed 
a rather constant f r . at low effective temperatures (see also 
Fig. 6 in iGuver et al]|2010aT) . as well as the fact that most of 
the short PRE bursts have a constant norm alization at late 
phases . We note here that the factors f c in iMaiczvna et alj 
(2005) correspond to the ratio of the energy where the peak 
of the model flux Fe is reached to the peak energy of the 
blackbody spectrum at effective temperature. Moreover, the 
low-luminosity models were calculated for high surface grav- 
ity instead of low effective temperatures, w hich leads to in- 
correct results (see Sulei manov et al.ll201 ll for details). The 
color-corrections obtained in this way, however, should not be 
compared to the data at all, because the color temperatures of 
the time-resolved spectra from X-ray bursts are computed by 
fitting the actual data in a specific energy interval (e.g. 3-20 
keV for RXTE/PCA) with the diluted blackbody function with 
arbitrary normalization. The values of f c shown in Fig. |6]on 
the other hand are produced using the procedure similar to 
that applied to the data, i.e. by fittin g the model spectra in the 
3-20 keV range. As was shown by Suleima nov et alJ (1201 ll) 
the color-correction has a flat part at I ~ 0.2-0.5 for most 
of the chemical compositions, but the actual value of f c de- 
pends on the hydrogen fraction, and, for example, for X = 1 
it is closer to 1.5 than to 1.4 (see Fig. [6j. At lower luminosi- 
ties f c can first drop because of iron edges and then increase 
to rather high values. Thus, there is no unique value for f c 
in the cooling tail. The expected significant variations of f c 
with flux also imply that constancy of the apparent blackbody 
area in the cooling tail contradicts the burst atmosphere mod- 
els and therefore the bursts showing such behaviour obviously 
demonstrate influence of some physics not included in these 
simplest models of NS atmospheres, and thus cannot be used 
for determination of NS masses and radii with the help of the 
aforementioned models. 

Third, different PRE bursts from the same source show dif- 
ferent cooling tracks, for example, the long burst and the short 
bursts of 4U 1724-307 (see Fig.HJl have normalizations dif- 
ferent by a factor of two. This makes the determination of the 
apparent area from a single burst not unique. 

And finally, the most serious problem with this approach 
is that out of the whole amount of information on the cool- 
ing tail of the burst, one uses only two numbers and it is not 
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Figure 8. Illustration of the suggested new cooling tail method. The depen- 
dence K~ 1 / A -F as observed during the cooling track of the long burst from 
4U 1724-307 on November 8, 1996 (circles). The theoretical f c -L/L Edd 
dependence is shown by the dashed curve (right and upper axes) and the 
best-fit relation (solid curve). 

checked whether these quantities are actually consistent with 
each other. For example, the theory also predicts that the color 
correction f c changes from wl.7 to ml A when the luminos- 
ity drops from the Eddington to about 1/3 of the peak value. 
This also implies that the blackbody normalization between 
the touchdown point and the decay phase must increase at 
least by a factor of two. It is really true for the long burst from 
4U 1724-307, while the short bursts have nearly constant K, 
which is two times smaller than that in the long burst, imply- 
ing probably a partial eclipse of the NS by the optically thick 
accretion disk and/or the influence of the boundary laye r on 
the structure of the NS atmosphere as discuss ed in Secti on|2JJ 
We also note here that al l bursts analyzed bv lOzel et al.l (2009) 
and lGiiver et al.l d2010allbh are short, they do not show enough 
variations of K in their cooling tracks, and therefore the re- 
sults obtained from these bursts are not reliable (see Sect. 15.21 
for more details). 

On the base of all these arguments we offer a new approach 
to the NS mass and radius estimations using the information 
from the whole cooling track. 

3.4. Determining M and R using the cooling tail method 

If the radiating surface area does not change during the 
burst decay phase, the evolution of the normalization is fully 
determined by the color correction variations (see Equa- 
tions (fTTT i and $Y2\ ). We thus suggest to fit the observed re- 
lation K^l^-F at the cooling phase of the burst by the the- 
oretical relations f r - L/L-f.m (shown in Fig. [6] see also 
ISuleimanov et al.ll201 lb with free parameters being A and the 
Eddington flux F^dd (see Fig. [8] for illustration). The be- 
haviour of / c depends rather weakly on the NS gravity and 
chemical composition, which substantially reduces the model 
dependence of the fitting procedure. Using the obtained best- 
fit parameters, we c an th en apply the method identical to that 
described in Section [3~3l 

The main advantages of the proposed cooling tail method 
is that there is no freedom in choosing f c in the cooling tail, 
the determination of the Eddington flux becomes decoupled 
from the uncertainties related to the touchdown flux as the 
whole cooling tail is used, and finally, one can immediately 
check whether the burst spectral evolution is consistent with 
theoretical models and whether the employed model includes 
the majority of the relevant physics for the description of the 
considered phenomenon. This check can help to choose for 
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Figure 9. Comparison of the X-ray burst data for 4U 1724-307 to the the- 
oretical models of NS atmosphere. The crosses present the observed depen- 
dence of K~ 1 ' 4 vs. F for the long burst, while diamonds represent two 
short bursts for the blackbody model with constant absorption /Yjj = 10 22 
cm" 2 . The solid curves correspond to the three best-fit theoretical models 
of various chemical compositions (se e Fi g. [6). The best-fit parameters i^dd 
and A, defined by Equations {6) and 1121 . are given in TablefT] 

further analysis only those bursts that follow the theory. 

4. RESULTS 
4.1. The long burst from 4U 1724-307 
4.1.1. Determining NS parameters using the cooling tail method 



Let us apply the method described in Section [3T4| for deter- 
mining NS mass and radius from the data on the long burst 
from 4U 1724-307. We fit the dependence of the normaliza- 
tion constant K on the observed flux F for the long burst by 
the theoretical curves / c - L/L-^dd computed for three chem- 
ical compositions. They give a good description of the data at 
intermediate fluxes for the data points to the right of the verti- 
cal dashed line (Fig. [9j, but below the touchdown, which we 
use for fitting. Close to the touchdown, significant deviations 
are probably caused by deviations from the plane-parallel at- 
mosphere and effects of the wind (thus the models are not re- 
liable). Strong deviations are also visible at low fluxes where 
the burst spectrum is probably modified by accretion. As- 
suming that i*Edd is actually reached at the touchdown con- 
tradicts the following evolution of the parameters during the 
cooling phase. The fits are better for the hydrogen-rich atmo- 
spheres. The results of the fitting for all considered chemical 
compositions of the NS atmosphere are presented in TableQ] 
The uncertainties in A and i*Edd are obtained with a bootstrap 
method. 

T aking the distance in the range 5.3-7.7 kpc (see Section 
12.21 ). we convert a distribution of -Ffedd and A using Monte- 
Carlo simulations to the distribution of M and R (Fig. [TOland 
TableQ]). The resulting contours are elongated, because of the 
uncertainty in distance, along the curves of constant Edding- 
ton temperature. The pure helium model atmospheres give a 
mass which is too small from the stellar evolution point of 
view. It is also below the mass-shedding limit if the star is ro- 
tating faster than at about 500 Hz. Pure hydrogen atmosphere 
models are consistent with the data only for D < 6 kpc, while 
for the atmosphere of solar composition the upper limit is 
7 kpc. The hydrogen rich atmosphere models give a lower 
limit on the stellar radius of 14 km independently of the metal 
abundance (see TableQ} for NS masses less than 2.3M Q , and 
smaller radii are allowed only for high NS masses. For the 
helium atmosphere, the solution shifts towards higher masses 
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Table 1 

Best Fit Parameters 



Atmosphere model ^Edd -A Tkdd.oo x 2 /d.o.f. M R 

(10" 7 ergs" 1 cm" 2 ) (km/10 kpc)- 1 /2 (10 7 K) (M Q ) (km) 

Hydrogen 0.525±0.025 0.170±0.001 1.64±0.02 5.0/5 1.9±0.4 / 2.45±0.15 14.7±0.8 / 11.7±1.3 

1.4 (fixed) 14.2±0.4 

Solar H/He, Z = 0.3Z© 0.521±0.020 0.172±0.002 1.66±0.02 5.8/5 1.85±0.6 / 2.7±0.15 15.5±1.5 / 13.0±1.0 

1.4 (fixed) 15.2±0.4 

Helium 0.50±0.02 0.178±0.002 1.71±0.02 11.3/5 l-05±^ 5 18.0±|| 

1.4 (fixed) 20.2±0.5 

Note. — Results of the fits to the K -1 ^ 4 -F dependence with the NS atmosphere models for various chemical compositions and log g — 14.0. For hydrogen 
and solar composition atmospheres there are two solutions for M and R (see Fig. 1101 . Neutron star mass and radius are computed from A and i*Edd assuming a flat 
distribution of the distance between 5.3 and 7.7 kpc with Gaussian tails of lcr=0.6 kpc. Errors correspond to the 90% confidence level. 




Neutron star radius R [km] 

Figure 10. Constraints on the mass and radius of the NS in 4U 1724—307 
from the long burst spectra (fitted with the blackbody model and constant 
absorption). The dotted curves correspond to the best-fit parameter A for 
the distance to the source of 5.3 kpc. For a flat distribution of the distance 
between 5.3 and 7.7 kpc with Gaussian tails of 1<t=0.6 kpc, the constraints 
are shown by contours (90% confidence level). They are elongated along 
the (dashed) curves cor responding to the Eddington temperatures TEdd,oo 
given by Equation U3> (which do not depend on the distance). These cor- 
respond to the three chemical compositions: green for pure hydrogen, blue 
for the solar ratio of H/He and subsolar metal abundance Z = 0.3Zq ap- 
propriate for Terzan 2 iOrtolani et al. 1997), and red for pure helium. The 
mass-radius relations for several equations of state of neutron and strange 
stars matter are shown by solid pink curves. The upper-left region is ex- 
cluded by constraints from the causality requirements lHaensel et alj |2007t 
Lattimer & Prakash 2007). The brown solid curves in the lower-right region 
correspond to the mass-shedding limit and delineate the zone forbidden for 
4U 1724-307, if it had a rotati onal frequency of 500 or 619 Hz, the highest 
detected for the X-ray bursters (Strohmaver & Bildsten 2006). 

and larger radii exceeding the mass-shedding limit (e.g. for a 
1.5 solar mass star the radius is about 20 km). If we take the 
canonical neutron star mass of 1.4M©, the NS radius is then 
strongly constrained at R = 15.2 ± 0.4 km assuming solar 
H/He composition with Z = 0.3Z Q and 14.2±0.4 km for hy- 
drogen. The obtained constraints (see Table [TJ imply a stiff 
equation of state of the NS matter. 

The choice of the highest-flux point at the K~ x l^-F plot 
used in the fitting procedure affects the results slightly. Ne- 
glecting even the second point after the touchdown, reduces 
the estimated FEdd by 3%, while A (defined by the horizon- 
tal part of the cooling track) remains unchanged. This leads 
to a 2% increase in radius and 4% decrease in the NS mass 
estimates. The fits to the K^^-F dependence, which is ob- 
tained with blackbody fits with free A^h, in the same flux inter- 
vals as for constant TVh give values of A about 10% smaller, 



while the Eddington flux estimates remain the same within 
1%. In this case, the estimated NS radii grow by 20% relative 
to those shown in Fig.fTol 

In addition, there is a systematic uncertainty of about 
10% in the absolute values of fluxes me asured by the RXTE 
dKirsch et al.ll2005l:IWeisskopf et al.l2010l) . It acts similarly to 
an additional 5% uncertainty on the distance and does not af- 
fect the value of TEdd,oo- With the current uncertainty on the 
distance to 4U 1724-307, this additional inaccuracy does not 
substantially increase the error bars on M and R in Table [T] 

The determined Eddington flux is smaller than the touch- 
down flux by about 15%. The main part of this differ- 
ence can be easily explained by the tempera ture depen- 
dence of the electron scattering opacity dPaczynskil 119831: 
iPozdniakov et al.lll983tPwin et al.l l 19931). In our model cal- 
culations (Sule imanov et al.ll201 ll) . we used the Kompaneets 
equation to describe the electron scattering and this assumes 
the Thomson opacity. As the upper atmosphere layers can be 
as hot as 3-3.5 keV, the electron scattering opacity there is 
smaller than the Thomson one by about 6-8% and the actual 
Eddington limit is reached at correspondently higher luminos- 
ity than iEdd given by Equation (0. Because we fit the data at 
luminosities much below the Eddington (where correction to 
the Thomson opacity are small), we determine FEdd as given 
by Equation ©, which is used for determination of M and R. 
We also note here that the opacity effect is not included in the 
touchdown method, which assumes the touchdown flux being 
equal to the Eddington flux at Thomson opacity. 

The remaining ^8% difference between the "true" Edding- 
ton flux and the touchdown flux can be understood, if we take 
into account that the maximum luminosity for PRE bursts can 
exceed the Eddington luminosity on the surface due to the 
dependence of the observed Edding ton lum inosity on the red- 
shift z (see Equation ([5]) and Lewi n et al.ll 1993b . This would 
imply that the photospheric radius corresponding to the touch- 
down exceeds the NS radius by ^25%. The corresponding 
color correction then has to be ~ 15% larger than for our mod- 
els with I = 0.98, which is co nsistent with that expected at 
L ~ L Ed d (lPavlovetalJU991h . The sharp maximum in Zbb 
(and minimum in K) can arise from a joint influence of the 
increasing color correction and decreasing effective tempera- 
ture during the photospheric radius expansion phases. 

The lower limit on the NS radius of 13.5-14 km as ob- 
tained by us is consistent with the measurements for the ther- 
mally emittin g quiescent NS X7 i n the globular cluster 47 Tuc 
(14.5ti g km: Hei nke et alJ2006T) . However, the radii of other 
thermal ly emitting quiescent NSs are significantly smaller (9- 
13 km; IWebb & Barrel 120071: iGuillot et alJ feOl 1). We note 
that these results depend on the model of the NS hydrogen at- 
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Figure 11. Comparison of the solutions (thick curves) obtained for the long 
burst using the cooling tail method (solid curves) and the touchdown method 
(dashed curves), and for the short bursts using the touchdown method only 
(dotted curves). The curves corresponding to the obtained TEdd.oo ar >d A (at 
-DiO = 0.47) are also shown. All solutions are derived for X =0.7374. The 
dash-dotted curve corresponds to the short burst with K = 115, i^dd -7 = 
0.55 and f c = 1.78 (i.e. T Edd>00 = 1.45 keV). 

mosphere , which were also c omputed for the passively cool- 
ing NSs (Hein ke et al.l 12006). In case a quasi-spherical ac- 
cretion occurs at a low rate during the quiescence, the addi- 
tional heat dissipated in the upper atmosphere wo uld increase 
the te mperature there, and, therefore, / c (see e.g. IZane et aT] 
2000). In this case, the NS radius can be underestimated. We 
also need to note that if the NS in 4U 1724-307 has a spin 
of 500 Hz, the radius of the non-rotating star would be about 
1 km smaller than estimated above. Our constraints on the 
NS radii are in good agreement with the NS radii (13-16 km) 
evaluated by Suleimanov & Poutanen (2006) from the spectra 
of low-mass X-ray binaries using the spreading layer model. 
Rather l arge NS radii are a llowed by the modern equations 
of state (Hebeler et al. 2010), which predict the upper limit of 
13.5 km for a 1.4M Q NS. 

4.1.2. Touchdown method 

We ap ply now the touchdown method described in Sec- 
tion [33] to the same data on the long burst. This approach is 
used for illustrative purpose only and we do not consider here 
any statistical errors. Let us take the Eddington flux equal to 
the touchdown flux .FEdd, -7=0. 605, the normalization in the 
cooling tail K=230, / c =1.4 and X=0.7374. From these quan- 
tities we obtain TEdd.oo = 1-84 x 10 7 K and an upper limit 
on the distance Dio,max = 0.5. The curves relevant to the de- 
rived TEdd,oo and the distance Dio = 0.47 (corresponding to 
la deviation from the minimal distance of 5.3 kpc) are shown 
in Fig.QT|by dashed lines. Using this approach we can esti- 
mate M (1.6 ± 0.2)M Q and R ss 10.0 ± 1.5 km. We see 
that the touchdown method gives a substantially smaller NS 
radius compared to the cooling tail method. 

The uncertainties in M and R are actually much larger, 
because of the uncertainties in chemical composition, color- 
correction and distance (see Fig.QJJ. The most significant 
errors arise due to unknown distance and the hydrogen mass 
fraction in the NS atmosphere (see top panel of Fig. [12] where 
we considered two limiting cases for the distance of Dm = 
0.47 and 0.83 corresponding to la deviations on both ends 
of the distance distribution and for the hydrogen mass frac- 
tion X = 1 and 0). Uncertainties in f c are also affecting 
the results (middle panel of Fig.[T2l. for example, changing 
f c from 1.35 to 1.45 increases the maximum possible R from 
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Figure 12. Expected uncertainties for the M and R solutions obtained for 
the long burst using the touchdown method. 

10 to 13 km and M from 1.5 to 2M©. We note that f c is 
actually closer to 1 .5 for the hydrogen-rich atmospheres (see 
Fig. [6}. Less significant errors appear due to uncertainties in 
the observed Eddington flux on the NS surface (bottom panel 
of Fig.[l2]i. A larger flux corresponds to smaller M and R. 

4.2. Short bursts 

The cooling tracks for the short PRE bursts shown in Fig. 
[9] are very different from that of the long burst and are com- 
pletely inconsistent with the theoretical dependences. This is 
a strong argument that these cooling tracks are affected by 
some additional physics and cannot be used for determination 
of NS mass and radius. Ignoring that fact, let us still apply 
the touchdown method to the short bursts. Taking from the 
observed data (see Fig. [4} FEdd, -7 ~ 0.55 and K « 115 
and assuming f c = 1.4, we find that the curves on the M—R 
plane corresponding separately to constraints from -FEdd and 
K do not cross for hydrogen-rich atmospheres (X > 0.7374) 
at any distance larger than 3.8 kpc, while for D — 4.7 kpc 
they cross only for X < 0.45 (see Figs.[TT]and[9]). The situ- 
ation that the two observables are consistent with each other 
in a very restricted range of X and distances is not unique 
to 4U 1724-307, bu t typical for m any bursters, par t icularly 
those analyzed b ylOzel et al.l d2009l) and lGiiver et alj d2010al) 
as was shown bv lSteiner et alj (120101) . see Section [5721 for de- 
tails. For pure He atmospheres, there are consistent solutions 
at M ~ 1.5M Q and R ~ 10 km at the distance ~ 6 kpc. 
However, we note again that the results from these short bursts 
are not reliable, because their cooling tracks contradict the 
theory the analysis is based on (which predicts f c ~ 1.4 in 
the cooling tail and much higher at luminosities close to the 
Eddington). 

The measured K ss 115 from the horizontal part of Fig. [4] 
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bottom panel is two times smaller than corresponding values 
for the long burst K « 230. This can result from two ef- 
fects simultaneously. The long X-ray burst occurred in the 
hard state, when the accretion flow in the NS vicinity had a 
small optical depth which only marginally affects the NS pho- 
tosphere and it cannot eclipse the neutron star. On the other 
hand, the short bursts happened in the soft state, when the per- 
sistent emission originated in the optically thick accretion disk 
and the boundary layer. Therefore, the accretion disk can just 
block part of the star in the decay phase of the burst reducing 
K by factor of up to two for the case of large inclination. Even 
if the inclination is small, the apparent NS area decreases by 
a factor of k = 1 — u 2 (where u — Rs/R). These limiting 
cases can be united in a simple approximate formula for the 
reduction factor 

k = - (l + [l - 2u 2 ] cosi) . (21) 

An additional effect can be related to the optically thick 
boundary layer. If the s preading layer model describ es the 
boundary layers correctly (Inogamov & Sunvaev 1999), a sig- 
nificant part of the emergent radiation can arise in the rapidly 
rotating spreading layer above the hot NS surface, which has 
a reduced effective gravity due to the centrifugal force, re- 
sulting in a flux through the atmosphere close to the local 
Eddin gton limit and a high color correction f c rj 1.6-1.8 
(see [Suleimanov & Poutanen 2006) and therefore small K. 
Using X -1 / 4 sa 0.305 as measured in the short bursts and 
A « 0.172 as determined from the long burst, we can estimate 
then the color-correction at the cooling tails of the short bursts 
f c sa 1.77 k 1 / 4 , which accounts for both effects (see Fig.fTTb. 
The importance of the spreading layer can vary with the ac- 
cretion rate and potentially the normalization val ues from K 
for th e long burst to K/2 can be found (see e.g. Zhang et al. 
l2OT0l) . 

It is also possible that the chemical composition of the NS 
atmosphere during short and long bursts is different. For ex- 
ample, pure helium atmospheres give K by 20-25% larger 
than pure hydrogen atmospheres because f c for helium at- 
mospheres is 5-6 % smaller (see Equation (fTTTi and Fig. [6}. 
However, the maximum temperature at the touchdown point 
must be « 20% larger for helium atmospheres (« 3.2 keV, 
see Equation (fT3l), while the maximum temperatures in the 
long and the short bursts are very close to each other, « 2.7 
- 2.8 keV, which argues against a difference in the chemical 
composition. 

5. COMPARISON TO OTHER X-RAY BURSTERS 
5.1. Long PRE bursts 

iPenninx et al.l d 19891) analyzed two long-duration (>100 s) 
PRE bursts observed by EXOSAT/ME in 1984 and 1986 from 
4U 1608-52 in its hard state at a rather low persistent flux of 
(l-2)xl0~ 9 erg cm -2 s" 1 in the 2-20 keV band. The evo- 
lution of Tbb-F 1-1 / 4 (which is proportional to K^ 1 / 4 ) with F 
shown in their Fi g. 7 is almost identical t o the models pre- 
s ented in Fig. 6 of Suleimano v et al.l (12011b . 

IKuulkers etafl d2003l) reported observations by RXTE/PCA 
of the long PRE bursts from 4U 1724-307 (analyzed here), 
the atoll source 4U2129+11 in globular cluster M15, and 
H1825-331 in globular cluster NGC6652. Spectral evo- 
lution after th e touch down in all sources is very similar. 
IKuulkers et al.l d2003l) also reported observations by Bep- 
poSAXTWFC of 24 long PRE bursts from4U 1724-307. Spec- 
tral evolution during the cooling phases of these bursts is 



entirely consistent with that observed in the long burst by 
RXTE/PCA. 

Observations by Ginga/LAC of a long PRE burst from 
4U 2129+11 during its island (hard) state at a low persis- 
ten t flux level of ~ . 5 x 1 ~ 9 erg cm -2 s _1 are presented 
by Ivan Paradijs etai] ( fl990h . The behaviour of T bb .F~ 1/4 
at fluxes above 30% of the peak (touchdown) flux shown in 
their Fig. 10 is very similar to that for the long burst from 
4U 1724-307 shown in our Fig. [9] For both objects the data 
at high fluxes are well described by the theory. We note that in 
both cases the position of the touchdown point in the K^ 1 / 4 - 
F diagram is not consistent with extrapolation of the data 
from intermediate fluxes, implying that the Eddington flux is 
smaller than the touchdown flux. 

5.2. Short PRE bursts from EXO 1745-248, 4U 1820-30 and 
4U 1608-52 

Recently strong claims appeared in the literature that both 
the NS mass and radius can be determined with an accu- 
racy of better than 10% from the PRE bursts from three NSs: 
EXO 1745-248 (M = 1.7±O.1M , R = 9±1 km,|OzeIetaD 
I2009I). 4U 1608-52 (M = 1.74 ± 0.14M©, R = 9.3 ± 1.0 
km. iGiiver etail 120103) . a nd 4U 1820-30 (M = 1.58 ± 
Q.O6M .i? = 9.1±0.41an. lGuveretaH 2010bl). The authors 
of these papers used only short PRE bursts, which, as we have 
seen above, are suspicious, because their spectral evolution is 
not consistent with the theory the method is based on, with 
the main reasons being probably partial blocking of the NS 
by the accretion disk and the effects of the spreading layer on 
the NS atmosphere. A high declared accuracy cannot be un- 
derstood in the light of all the uncertainties, especially on the 
distance and the chemical composition (see Fig. [121). We try to 
find below the answers by the critical consideration of the in- 
put numbers and the assumptions made in the aforementioned 
papers. 

5.2.1. EXO 1745-248 in Terzan 5 

lOzel et ail d2009l) have determined the following param- 
eters from two PRE bursts from EXO 1745-248: F Ed d.-7 
= 0.625±0.02 and K = 116±13. These correspond to 
Tkdd.oo = 2.2 x 10 7 K and the maximum possible distance 
(see Equation d20l i) -Dio.max = 0.5978 at / c =1.4 (which was 
fixed). For the chemical composition, the authors also as- 
sumed pure helium, X = 0. The distance was taken Dio = 
0.63±0.0315 (box-car distribution) with the strict lower limit 
of 0.5985 being very close (within 0.1%) to the maximum 
possible distance for the observables. As a result the curves 
corresponding to TEdd.oo and K only touch each other (see 
Fig-US. We note that the assumption of the box-car distribu- 
tion for the distance (as well as for K) only allow s distances 
within 10% above the used strict lower limit. Thus lOzel et alJ 
(2009) de facto assumed nearly the delta-function distribution 
of the distance at ^6 kpc. Fixing also f c and X they thus com- 
pletely removed all uncertainties connected with these param- 
eters. The declared errors in M and R thus reflect the statis- 
tical errors in i*Edd,-7 and K only, which are of course small 
bec ause of brightness of the considered events. We also note 
that lOrtolani et al j (120071) evaluated the distance to Terzan 5 of 
5.5±0.9 kpc and lValenti et all (120071) gives D = 5.9 kpc. Re- 
laxing the assumption of the box-car distribution for the dis- 
tance (e.g. by using a Gaussian distribution) will inevitably 
move the solutions towards smaller distances and therefore 
smaller masses and radii. 
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Figure 13. Curves corresponding to constant TEdd.oo ar >d K (for the strict 
lower limit on the distance) at / c =1.4 and X = for EXO 1745-248 (thick 
solid curves), 4U 1820-30 (thick dashed curves), and 4U 1608-52 (thick dot- 
ted curves). The contours for constant Tgdd oo at /c=1.5 and X = 0.7374 
for double normalization IK are also presented (thin curves). 

5.2.2. 4U 1820-30 in NGC 6624 

Using the RXTE data f or five short PRE bursts from 
4U 1820-30, iGtiver et al.l d2010bl) determined F Ed d,-7 = 
0.539±0.012and.rY = 92±2. These correspond to T Edd oo = 
2.25 x 10 7 K and the maximum possible distance -Dio,max = 
0.617 at f c = 1.4 and X = 0. These authors fixed the chem- 
ical composition at X = and varied f c between 1.3 and 
1 .4 with the strict limit s (box -car distribution). Using esti- 
mates bv lKuulkers et al.l (2003) and lValenti et all (120071) . they 
took the distance to the source Diq = 0.82±0.14, but again 
assumed strict limits (box-car distribution). 

Formally, there is no solution possible for these observ- 
ables as the curves corresponding to TEdd.oo and K do not 
cross (see Fig.[T3~l>. The solutions for the larger distances in 
the Gaussian tails of the distributions of i^dd -7 and K are 
still p ossible with the probability of about 10 _,> (ISteiner et al.l 
2010). Again the obtained errors in M and R reflect the statis- 
tical errors in -FEdd,-7 and K only, but the solution is highly 
unlikely. 

As in the case of EXO 1745-248, relaxing the distance con- 
straints to allow smaller distances moves the solution towards 
smaller distances and smaller M and R. Similarly to the case 
of the short bursts from 4U 1724-307, the spectral evolution 
of the considered bursts is not consistent with theory. The 
value of K^ 1 / 4 drops by only 12% after touchdown, while 
the theory predicts 20% variations. 

5.2.3. 4U 1608-52 

Using the data fr om four PRE bursts a nd one non-PRE burst 
from 4U 1608-52, IGtiver et al.1 d2010al) found the following 
observables: F Ed d,-7 = 1.541±0.065 and K = 324.6±2.4. 
These correspond to TEdd,oo = 2.13 x 10 7 K and the maxi- 
mum possible distance £*io, ma x = 0.405 at / c =1.4 and X=Q. 
The chemical composition was allowed to vary from X = 
to 0.7, and / c within the interval between 1.3 and 1.4 with the 
strict limits. The distance was taken by these authors as Dio 
= 0.58^° ^ 8 with the strict lower limit equal to 3.9 kpc. This 
lower limit in this case is smaller than the maximum possible 
distance and the solutions exist (the curves corresponding to 
Tkdd.oo and K do cross, see Fig.[T3l. Formally, the authors 
allow variation of X in a wide range, but in reality only solu- 
tions with X < 0.04 are possible. At larger X the maximum 
possible distance becomes smaller than the strict lower limit 



on the distance of 3.9 kpc. As in the previous two cases, re- 
laxing the distance constraints to allow for smaller distances 
will move the solution towards smaller distances and smaller 
M and R. 

We a lso note here that bursts selected by IGtiver et"afl 
(2010a) occurred at a high accretion rate in the soft state. The 
normalization K is nearly constant during the cooling tails, 
which strongly contradicts the theory. On the other hand, the 
PRE bursts happening at a low accretion rate follow spectral 
evolution predicted by the theory, have normalization K about 
60% larger, and, of course, the NS mass and radius deter- 
mined from these data are also different (Poutanen et al., in 
preparation). 

5.2.4. Summary 

We conclude that small uncertainties in M and R obtained 
for the three bursters are the direct consequence of fixing the 
color correction / c , hydrogen mass fraction X, and most im- 
portantly of the assumed strict lower limit in the distance dis- 
tribution, which allows only solutions at its lower edge. The 
lower limits on the distance assumed by the authors in two 
cases turned out to be very close to the maximum possible dis- 
tance allowed by the observables. This completely removed 
all the uncertainties in the distance. Assuming a Gaussian 
distribution in distance allows solutions with much smaller 
(probably unphysical) M and R for all three sources. 

The results of the spectral fitting for the short PRE bursts 
for these three NSs are very close to the results obtained for 
the short bursts of 4U 1724-307. We suggest that during short 
bursts roughly half of the NS is visible most probably due to 
the eclipse by the accretion disk. Taking the apparent area 
twice the observed one, the determined NS radii move to the 
values consistent with those determined from the long burst 
from 4U 1724-307 (see Fig. Il3l . Additional corrections are 
also possibly needed because of the influence of the bound- 
ary layer on the dynamics and the spectral properties of the 
bursting atmosphere, resulting in a higher color correction 
f c ~ 1.6-1.8. 

Recently Stei ner et al.1 feOlO) suggested that the absence of 
the solutions for short bursts mentioned above can be fixed 
by relaxing the assumption that the Eddington flux is reached 
at the moment of touchdown (or rather that the photosphere 
at touchdown is not at the NS surface). This, however, does 
not solve the problem because the determined blackbody nor- 
malization in the cooling tail and the Eddington flux are not 
reliable in these short PRE bursts and the assumed distance 
constraints were too strict. As the spectral evolution during 
the cooling tails of the short PRE bursts (at high persistent lu- 
minosity) strongly contradicts the theory, the results obtained 
from these data are questionable. For the determination of the 
NS masses and radii we suggest to use only those data that do 
follow the theory. 

6. CONCLUSIONS 

We applied a recently computed detailed set of NS at- 
mosphere models cover ing a large range of luminosities 
dSuleimanov et al.l |20TT1) to the data of the PRE bursts of 
4U 1724-307. We showed that the variation of the apparent 
blackbody radius during the cooling stage of the 150 s long 
PRE X-ray bursts in 4U 1724-307 is entirely consistent with 
the theoretical color-correction-flux dependence at interme- 
diate fluxes. We thus obtained the Eddington flux and the ap- 
parent NS radius (divided by the distance to the source). We 
find that the Eddington flux (for Thomson opacity) is reached 
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not at the so called "touchdown", but later at a 15% lower 
luminosity. We constrained mass and radius of the NS using 
the estimated distance to the source. We find a lower limit 
on the stellar radius of —14 km for M < 2.3M at 90% 
confidence independently of chemical composition. (If the 
NS in 4U 1724-307 has a spin of 500 Hz, the radius of the 
non-rotating star would be about 1 km smaller.) Smaller radii 
are possible only for more massive NS. These results support 
a stiff equation of state of the NS matter. We showed that 
hydrogen rich accreting matter is preferred and obtained an 
upper limit on the distance to 4U 1724-307 of about 7 kpc. 

We have also demonstrated that the cooling tracks of the 
short PRE bursts from 4U 1724-307 that occurred during the 
high/soft state do not follow the evolution expected from the- 
ory and the NS apparent areas are a factor of two smaller 
than that for the long burst. The probable reason is the par- 
tial eclipse of the NS surface by the optically thick accretion 
disk. An additional spectral hardening during the cooling tails 
associated with the influence of the boundary/spreading layer 
can be also important. Therefore, the constraints on the NS 
mass and radius obtained from such bursts are not reliable. 

Finally, we showed that previous analyses of the short PRE 
bursts from three sources EXO 1745-248, 4U 1820-30 and 
4U 1608-52 are questionable, because they ignore the fact 
that spectral evolution during the bursts is not consistent with 
the theory for the passively cooling unobscured NS that is the 
base for the analysis. Assuming that the touchdown flux is 
reac hed when the photo sphere is detached from the NS sur- 
face (ISteiner et al.ll2010h does not solve the problem because 
the determined blackbody normalization in the cooling tail 
and the Eddington flux are not reliable in these short bursts. 

We suggest that only PRE bursts showing spectral evolu- 
tion consistent with the theory should be used when estimat- 
ing NS masses and radii. Further improvement in accuracy of 
determination of the NS parameters requires new atmosphere 
models with the exact Compton scattering kernel, as well as 
accounting for the Doppler effect due to the rapid NS rota- 
tion. Rotation also introduces additional distortion to the NS 
spectra, because of the difference in the effective gravity at 
the stellar poles and the equator due to the centrifugal force. 
We plan to investigate the importance of these effects in future 
work. 
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